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Abstract 

We discuss the discrete spectrum of N particles in a curved planar wave- 
guide. If they are neutral fermions, the maximum number of particles which 
the waveguide can bind is given by a one-particle Birman-Schwinger bound in 
combination with the Pauli principle. On the other hand, if they are charged, 
e.g., electrons in a bent quantum wire, the Coulomb repulsion plays a crucial 
role. We prove a sufficient condition under which the discrete spectrum of 
such a system is empty. 

1 Introduction 

A rapid progress of mesoscopic physics brought, in particular, interesting new prob- 
lems concerning relations between geometry and spectral properties of quantum 
Hamiltonians. They involve models of quantum wires, dots, and similar systems. 
While in reality these are rather complicated systems composed of different semicon- 
ductor materials, experience tells us that their basic features can be explained using 
simple models in which electrons (regarded as free particles with an effective mass) 
are supposed to be confined to an appropriate spatial region, either by a potential 
or by a hard wall. A brief description of this approximation with a guide to further 
reading is given in Ref. fl]. In addition, such models apply not only to electrons in 
semiconductor microstructures; a different example is represented by atoms trapped 
in hollow optical fibers ||. 

It is natural that most theoretical results up to date refer to the case of a single 
particle in the confinement. On the other hand, from the practical point of view it 
is rather an exception than a rule that an experimentalist is able to isolate a single 
electron or atom, and therefore many-body problems in this setting are of interest. 
For instance, two-dimensional quantum dots which can be regarded as artificial 
atoms have been studied recently, usually in presence of a magnetic field, either 
for a pair of electrons or in the semiclasical situation when a Thomas-Fermi-type 
approach is applicable — cf. [3-6] and references therein. 
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In these studies, however, geometry of the dot played a little role, because the 
confinement was realized by a harmonic potential or a circular hard wall. This is not 
the case for open systems modelling quantum wires where a deformation of a straight 
channel is needed to produce nontrivial spectral properties. In particular, a quantum 
waveguide exhibit bound states if it is bent 0, 0, ||, protruded [9-11] or allowing a 
leak to another duct [12-14], and the discrete spectrum depends substantially on the 
shape of the channel. With few exceptions such as Ref. ||15|| , however, the known 
results refer to the one-particle case. 

It is the aim of the present paper to initiate a rigorous investigation of many- 
particle effects in quantum waveguides. We are going to discuss here a system of N 
particles in a bent planar Dirichlet tube, i.e., a hard-wall channel, and ask whether 
iV-particle bound states exist for a given geometry. After collecting the necessary 
preliminaries in the next section, we shall derive first in Section 3 a simple bound 
for the neutral case which follows from the Birman-Schwinger estimate of the one- 
particle Hamiltonian in combination with the Pauli principle. 

The main result of the paper is formulated and proved in Section 4. It concerns 
the physically interesting case of charged particles; the example we have in mind 
is, of course, electrons in a bent semiconductor quantum wire. The electrostatic 
repulsion makes spectral analysis of the corresponding Hamiltonian considerably 
more complicated. Using variational technique borrowed from atomic physics, we 
derive here a sufficient condition under which the discrete spectrum is empty. The 
condition is satisfied for N large enough and represents an implicit equation for the 
maximum number of charged particles which a waveguide of a given curvature and 
width can bind. Some other aspects of the result and open questions are discussed 
briefly in the concluding section. 



2 Preliminaries 

The waveguide in question will be modelled by a curved planar strip E in M 2 , of 
a constant width d = 2a . It can be obtained by transporting the perpendicular 
interval [—a, a] along the curve T which is the axis of £ . Up to Euclidean trans- 
formations, the strip is uniquely characterized by its halfwidth a and the (signed) 
curvature s h- > j(s) of T , where s denotes the arc length. We adopt the regularity 
assumptions of Refs. |1|, [/J: 

(i) Q is not self-intersecting, 

(a) alMloo < i , 

(Hi) 7 is piecewise C 2 with 7', 7" bounded, 

and restrict our attention to the case when the tube is curved in a bounded region 
only: 

(iv) there is b > such that 7(3) = for \s\ > b; without loss of generality we 
may assume that 2b > a . 
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As usual we put % = 2m = 1 ; then the one-particle Hamiltonian of such a wave 
guide is the Dirichlet Laplacian — defined in the conventional way — cf. |L6 



Sec. XIII. 15. Using the natural locally orthogonal curvilinear coordinates s,a in E 
one can map — A£ unitarily onto the operator 

H x = -d s (l + U1 )- 2 d s - d 2 u + V(s,u) (2.1) 

on L 2 (M x (—a, a) ) with the effective curvature-induced potential 

vis u) = ^ 2 + _ 5 «v(«) a (2 2) 

1 ' ; ' 4(1 + m 7 (s)) 2 2(1 + m 7 (s)) 3 4(l + u 7 (s)) 4 1 ' ; 

which is e.s.a. on the core D(H) — {ip : ip 6 C 00 , ip(s,±a) = 0, if?/> G L 2 } - 
c/. Refs. 0] for more details. 

If the waveguide contains iV particles, the state Hilbert space is L 2 (S)) 7V ; the 
Pauli principle will be taken into account later. We assume that each particle has 
the charge e ; using the same "straightening" transformation we are then able to 
rewrite the Hamiltonian as 

N 

H N = H N (j, a, e) = J2 { ( l + u n( s j) Y 2 d* s ~ d l 3 + v ( s j' u j) } 

+ e 2 fi-nl- 1 , (2-3) 

l<j<l<N 

with the domain (H 2 (1R) ® Hq(—cl, a)) , where r} = fj(sj,Uj) are the Cartesian 
coordinates of the N-ih particle. 

As we have said our main aim in this paper is to estimate the maximum number 
of particles which a curved waveguide with given 7, a can bind, i.e., to find condi- 
tions under which the discrete spectrum of Hn is empty. To this end, one has to 
determine first the bottom of the essential spectrum. In complete analogy with the 
usual HVZ theorem pi, we find 



C"css 

(H N ) - 

where /z/v-i := inf a(Hff-i) ■ Obviously, 



7T ^ 2 



),oo , (2.4) 



inf a css (H N ) < fi N _ k + k (^-\ 



holds for k — 1 N — l , so 



inf a css (H N ) < N (t^) 2 ■ ( 2 - 5 ) 

In a straight tube the two expressions equal each other, while for 7 ^ we have a 
sharp inequality because < \^) holds in this case. 
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3 Neutral fermions 



If the particles in question are neutral fermions, one can get a simple upper bound on 
the number of bound states using the one-particle Hamiltonian fl2.1| ); it is sufficient 
to estimate the dimension of (Tdisc(-ffi) and to employ the Pauli principle. To this 
aim, one has to estimate H\ from above by an operator with the transverse and 
longitudinal variables decoupled; its projections to transverse modes are then one- 
dimensional Schrodinger operators to which the modified Birman-Schwinger bound 
may be applied [17-19]. In Ref. H] we used this argument in the situation where a 
is small so that only the lowest transverse mode and the leading term in Q2.2Q may 
be taken into account. 

A modification to the more general case is straightforward. We introduce the 
function 

^ obl ^ 
W{S) A51 + 25 s + A5i ' {6 ] 

where 

5± := liajHIoo, (3.2) 
which majorizes the effective potential, V(s,u) < W(s) . Furthermore, we set 

Wfc) := max jo, Ql) 2 (!-/»)} (3.3) 

for j — 2, 3, . . . ; in view of the assumptions (%%), (Hi) only finite number of them is 
different from zero. 

Replacing K by IV, and (l + wy)~ 2 by <5+ 2 , we get an estimating operator 
with separating variables, or in other words, a family of shifted one-dimensional 
Schrodinger operators; we are looking for the number of their eigenvalues below 

inf <7 ess (.Hi) = (-^j . The mentioned modification of the Birman-Schwinger bound 
is based on splitting the rank-one operator corresponding to the singularity of the 
resolvent kernel £~ e _K l s_s I at k = and applying a Hilbert-Schmidt estimate to 
the rest. In analogy with Refs. f[7 , [T8| , |19| we employ this trick for the lowest-mode 



component of the estimating operator, while for the higher modes we use the full 
resolvent at the values Kj := y / j 2 — 1 . In this way we arrive at the following 
conclusion: 

Proposition 3.1 The number iV of neutral particles of half-integer spin S which 
a curved quantum waveguide can bind satisfies the inequality 

n < (Wi + ^kMtffl^ 

a6 2 



+ / W^ds}. (3.4) 
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Remarks 3.2 (a) As we have said, the number of nonzero term in the last sum is 

finite. More exactly, the index j runs up to the entire part of y^l + {?jfj Halloo ; 
hence if a is small enough this term is missing at all. 

(b) The assumption (iv) is not needed here. It is sufficient, e.g., that the functions 

7,7', and ^"l 1 ^ 2 decay as |s|~ 1-e as \s\ — > 00. 



4 Main result: N charged particles 

We have said in the introduction that the present study is motivated mainly by 
the need to describe electrons in curved quantum wires. Unfortunately, the above 
simple estimate have no straightforward consequences for the situation when the 
particles are charged. While the electrostatic repulsion adds a positive term to 
the Hamiltonian ( |2.3| ), it may move at the same time the bottom of the essential 
spectrum since the energies of the bound "clusters" are, of course, sensitive to the 
interaction change. 

We need therefore another approach which would allow to take the repulsion 
term in (|2.3|) into account. An inspiration can be found in analysis of atomic N- 
body Hamiltonians. To formulate the result we need some notation. Given a positive 
(3 we denote by {A m }^ =1 the ordered sequence of eigenvalues of Dirichlet Laplacian 
at the rectangle 

~P5 + , 1(35+ 



R 

and set 



2' 2' 



x [—a, a) , (4.1) 



2E^ =i A m ... N = 2n 

UN) := { (4.2) 
2E^=iA m + A n+1 ... N = 2n + 1 

We have in mind here electrons and assume that the spin is ~ , otherwise Tp(N) 
has to be replaced by the sum of the first iV eigenvalues of 2S+1 identical copies 
of the Laplacian. Now we are able state our main result: 

Theorem 4.1 Assume (i)-(iv). o"di SC (#jv(7, a , e )) = for > 2 if the condition 

T ° [N) + ^T7 mN - 1] ^ ll#IUJV + (£) N + WT2 <43) 
is valid for some (3 > max{ 2b, 596 e~ 2 } . 

Proof: We use a variational argument which relies on a suitable decomposition of 
the configuration space. Consider a pair of smooth functions v, g from M + to [0, 1] 
such that 

... * < 1 
v(t) = { (4.4) 

1 ... t>§ 
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and 

v{tf+g(tf= 1. (4.5) 

Elements of the configuration space are (s,u) with s = {si,...,sjv} and u = 
{ui, . . . , un} . We denote ||s||oo := niaxjsi, . . . , sjy} and employ the functions 

s >- vQlsWnP- 1 ), gihUp- 1 ) , 

where (3 > 2b > a is a parameter to be specified later. By abuse of notation, we use 
the symbols v, g again both for these functions and the corresponding operators of 
multiplication. It is straightforward to evaluate ([H^, v]tp, vip) and the analogous 
expression with v replaced by g for a vector ip e D(H^) ; in both cases it is only 
the longitudinal kinetic part in fl2.3|) which contributes. This yields the identity 



(H N ip,ip) = (H N vtp,vilj) + (H N g%l),gip) 



N { II 2 - 2 1 

7=1 J 



where we have used the shorthands Vj := , gj := J^- , and 7j := 7(sj) . Notice 

further that the factors (l+u^-) -1 may be neglected, because Vj gj are nonzero 
only if Sj > (3 > 2b in which case jj — . Furthermore, with the exception of the 
hyperplanes where two or more coordinates coincide (which is a zero measure set) 
the norm ||s||oo coincides with just one of the coordinates s\, . . . , s n , and therefore 

N 

£{lMf+ IMH 2 } < W 2 max {\\ Vj \\l+ \\g 3 \\l} < /r 3 C |M| 2 , (4.6) 
j=i - 3 - 

where C := ||f / ||^ + \\g'\\lo ■ We arrive at the estimate 

(H N ip,ip) > L!M+A[#] (4.7) 



with 



LM := {H N <PA) ~ ^11*,, (4-8) 



where the last index symbolizes the norm of the vector restricted to the subset 
Hp := | s : /? < || s || oo < ^ | of the configuration space. 

Next one has to estimate separately the contributions from the inner and outer 
parts. Let us begin with the exterior. We introduce the following functions: 



Ms) = ^(2 Sl || S ||- 1 ) , 

i-i 

fj(s) = ^ (2^-lkH- 1 ) II ^ (2^|| S ||- 1 ) , j=2,...,iV-l 

n=l 

iV-1 

fM = U 9 (2Sn\\sV] 



n=l 
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It is clear from the construction that 



N 



Moreover, the functions 



E/;(*) 2 = !■ 



oo / 



(4.9) 



have a non-zero derivative only if \sj\ > IIHI^ 1 - Hence on the support of s i— > 
'KINloo/^ -1 ) the derivative is non-zero if \sj\ > \(5 > b . In other words, the function 
s l— *■ /i( s ) 2,y (ll s lloo/3~ 1 ) has zero derivative in all the parts of the configuration space 
where at least one of the electrons dwells in the curved part of the waveguide. 
Commuting the (longitudinal kinetic part of) Hn with fj , we get in the same way 
as above the identity 



N 



iiW = E{ L i[/^]-ll(v s />vil 2 } , 



(4.10) 



where V s := (d Sl , . . . , d Sl ) . Next we need a pointwise upper bound on Z^=i(V s /j) 2 : 
denoting <Xj := 2s_ 7 -||s|| 00 , we can write 



El(vj,)(«)l 2 = 



+ ( ? / ^l)^(^) 2 + ^(^l)V( ( 7 2 ) 2 + --- 

+#Vi) 2 <?(^) 2 . • • g(a N ) 2 + ■■■ + g^) 2 . . . g(a N ^) 2 g\a N ) 2 j , 
which gives after a partial resummation 

( v '(a 1 ) 2 + g'(a 1 ) 2 + g(a 1 ) 2 g'(a 2 ) 2 + 



-g{o l ) 2 ...g{o N - 1 ) 2 g\o N ) 2 } 



< 



V 



l S 'loo I j = l I II s lie-. 



recall that Co := ||'i/|| 2 K) + , x 



'" 2 II "'II- 2 . Consequently, 

AT 



1-1 



E^i[/W>]-4ATCo /^Hall" 1 



E L 2[/j^] , 
J'=l 



(4.11) 
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where 

L 2 [0] := Li[0] - ^Col^lMI- 1 ! • (412) 
Hence we have to find a lower bound to L^i^j) with := fjVip) . Since Sj > 
||| s lloo > |/3 > b holds on the support of ipj , we have V(sj,Uj) = there. This 
allows us to write 

2 2 ^ / _ x 

where H^_i refers to the system with the j—th electron excluded, and therefore 



mi 2 + e 2 £ nrVi.V'j 
#1=1 



Since |f^ — rj| < J(sj — si) 2 + Aa 2 < 2 J\\s\\l + a 2 , we have 



V2a/ 



Ml 2 + 



I 2 -U „2W/2 



The sought lower bound then follows from ( 4.12| ) and ( 4.8 



L 2 [^j] > Vn-i 



- c p~ 



(,2a) 



ANC Q ^ 3 \\s\\-J 



e 2 {N-l) 



recall that j\Tp := { s : /3 < ||s||oo < if} ■ The second and the third term at the rhs 
can be combined using 

4ivco||^INI|- 1 || 2 + cyr 2 !!^ < (4iv+i)Co||^|| s ||~ 1 || 2 . 

Furthermore, ||s||oo > [3 > 2b > a yields (||s|| 2 + a 2 ) 1 / 2 < \/2||s||oo and 



L 2 [4>j} > fJ-N-l 



( 7T V 



+ 



\2a/ , 
e 2 (iV-l) C (4iV+l) 



2v/2 







1-1 

Oll"lloc 



(4.13) 



We are interested in the situation when the second term at the rhs is positive. This 
is achieved if 

e 2 (iV-l) C (4iV+i) 
2V2 
which is ensured if we choose (3 in such a way that 

18^0, 



P > 



(4.14) 



recall that N >2 . Owing to the identity ( [4.1 Ij ) we then have 



7T ^ 



£iM > ( HN-1+ (^) ] IM| 2 , (4.15) 



which means in view of ( |2.4j ) that the external part of ip does not contribute to the 
discrete spectrum. 

Let us turn now to the inner part. The corresponding quadratic form in the 
decomposition ( (4.7|) can be estimated with the help of ( |2.3|) and ( |4.8|) by 

> s- 2 \\v s9 i>\\ 2 + ||v„#|| 2 + 



l«k<N 



+ e 2 £ (||r,-rl||-V,#) - pll#f; (4.16) 



recall that <5 + := l + a^Hoo. Using the function W defined by Q3.1p we find 

in«i>%)i<^(^-),so 

max{U(s,w) : (s,u) e M x [—a,a]} < \\W\loo . 
Consequently, the curvature-induced potential term can be estimated by 

N 

E {v{ Sj , Uj )g^9ip) < niyiuivii^ii 2 . 

3=1 

Furthermore, on the support of g we have 



\rj-f k \ < 2 a 2 < pf3 2 +^a\ 

because ||s||oo < |/3 holds there. At the same time, (3 > 2b > a , so we arrive at 
the estimate 

\rj-r k \ < V7p, 

which yields 

2^ [\\ r j- r k\\ 9^,9^) > 



l«k<N 



Now we can combine the above estimates with the inequality ^ < j^y^ which 
follows from (|4. 14j ) to get the bound 



LiM > <5; 2 ||v s ^|| 2 + ||V U ^II 2 

e 2 N(N-l) e 2 



+ 



-N\\W\\ oc + 



2fJV7 l8pV2 



imp- (4-i7) 
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Now we can put the above results together. In view of the inequality (|4.15| ) and of 
(|2.5|), the last bound tells us that if at has no discrete spectrum for N > 2 provided 



Sl 2 \\V s g^ 



|V U #|| 2 + 



e 2 N(N-l) 



- iV||W||oo- N 
for some (3 which satisfies the condition 



7T 

2a 



18/3^2 



> 



(4.18) 



P > max { 26, 



18V2C 



(4.19) 



The first two terms in ( [4.18j ) are nothing else than the quadratic form of the 2N- 
dimensional Laplacian on — cf. ( f4.1|) . By Pauli principle each eigenvalue may 
appear only twice, thus one has to take the orthogonal sum of two copies of the 
Laplacian on Rp and to summ the first N eigenvalues of such an operator. This is 
exactly the quantity which we have called Tp(N) . 

To finish the proof, it remains to estimate Cq which appears in the conditions 
(|4.14j) and (f4.19|) . We will not attempt an optimal bound and put simply 



for t-1 



v{£) := sin(V 2 (l-2£ 2 

f e (0, 1) , then 
has the maximum value 2v / 2(87r) 2 /3 ~ 595.5 . 1 



5 Conclusions 

Since the present study is rather a foray into an unchartered territory, the result is 
naturally far from optimal. Let us add a few remarks. First of all, it is clear that the 
overall size of the curved region affects substantially the number of particles which 
the waveguide can bind. We know that any curved tube has a one-particle bound 
state |TJ, hence a tube with N slight bends which very far from each other (so 
far that the repulsion is much smaller that the gap between the bound state energy 
and the continuum) can certainly bind iV particles for iV arbitrarily large. 

The method we use is borrowed from atomic physics where it yields bounds on 
ionization of an atom. Of course, there are differences. The binding is due to the 
curved hard wall of the waveguide rather than by the electrostatic attraction to the 
nucleus, and the spectrum of our one-particle operator (|2.1|) is finite. Consequently, 
there is a maximum number of particles which a given curved tube can bind as long 
as the particles are fermions. Bosons can occupy naturally a single state, and the 
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idea of a Bose condensate of neutral spin-zero atoms in a curved hollow optical fiber 
is rather appealing. 

On the other hand, a non-zero particle charge changes the picture, and even the 
number of bosons bind by a curved tube is limited: notice that the condition ( [4. 3D 
is satisfied for large enough N without respect to the Pauli-principle term Tp(N) . 
Of course, the fermionic nature reduces the maximum number N further, since 
Tp(N) growth for large N is between o(N 3 ) in the limit a — > and o(N 2 ) for 
2b ~ a . At the same time, the maximum number also depends on the value of the 
charge. Since — > and the remaining terms in ( |4.3| ) are independent od 
e we see that o~disc(H N ) = for any N > 2 provided e is large enough. Thus our 
result confirms the natural expectation that for a given curved tube and sufficiently 
charged particles just one-particle bound states can survive. 

We have not addressed in this paper the question about the minimum number of 
particles which a curved quantum waveguide can bind. The gap between the trivial 
result which follows from the one-particle theory [7], |J and the condition ( [4.3| ) 
leaves a lot of space for improvements. Moreover, it is a natural question whether 
strongly curved tubes which can bind many particles allow for some semiclassical 
description analogous to the case of the quantum dots || . This is a task for a future 
work. 
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